Spontaneous fission modes and lifetimes of super-heavy elements in the nuclear 

density functional theory 
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Lifetimes of super-heavy (SH) nuclei are primarily governed by alpha decay and spontaneous 
fission (SF). Here we study the competing decay modes of even-even SH isotopes with 108 < 
Z < 126 and 148 < N < 188 using the state-of-the-art self-consistent nuclear density functional 
theory framework capable of describing the competition between nuclear attraction and electrostatic 
repulsion. The collective mass tensor of the fissioning superfluid nucleus is computed by means of the 
cranking approximation to the adiabatic time-dependent Hartree-Fock-Bogoliubov approach. Along 
the path to fission, our calculations allow for the simultaneous breaking of axial and space inversion 
symmetries; this may result in lowering SF lifetimes by more than seven orders of magnitude in 
some cases. We predict two competing SF modes: reflection-symmetric and reflection-asymmetric. 
The shortest- lived SH isotopes decay by SF; they are expected to lie in a narrow corridor formed by 
280 Hs, 284 F1, and if|Uuo that separates the regions of SH nuclei synthesized in "cold fusion" and 
"hot fusion" reactions. The region of long-lived SH nuclei is expected to be centered on 294 Ds with 
a total half- life of ~1.5 days. 

PACS numbers: 24.75.+i, 25.85.Ca, 21.60.Jz, 27.90.+b, 23.60.+e 



Introduction — The SH nuclei represent the limit of nu- 
clear mass and charge; they inhabit the remote corner 
of the nuclear landscape whose extent is presently un- 
known. The mere existence of long-lived SH isotopes has 
been a fundamental question in science since Seaborg and 
Swiatecki coined the notion of the SH "island of stability" 

0. 

Theoretically, it is anticipated that the majority of SH 
nuclei would fission and/or a-decay, but predictions vary 
from model to model, primarily due to our inability to 
make accurate predictions of SF half-lives. Here the main 
uncertainty is our imperfect knowledge of effective nu- 
clear interactions and the strong Coulomb frustration ef- 
fects due to the interplay between the long-ranged elec- 
trostatic repulsion and the short-ranged nuclear force. 
By the end of the 1960s, it had been concluded that 
the existence of the heaviest nuclei with Z > 104 was 
primarily determined by the quantum-mechanical shell 
effects (i.e., single-particle motion of protons and neu- 
trons in quantum orbits) 0, Q . These early microscopic- 
macroscopic (MM) calculations predicted the nucleus 
with Z = 114, N = 184 to be the centre of an island of 
long-lived SH nuclei. This result stayed practically un- 
challenged until the late 1990s when self-consistent mean- 
field (SMF) models, based on realistic effective interac- 
tions, were applied to SH nuclei [4]. Currently, most the- 
ories agree that nuclei around N = 184 and Z between 
114 and 126 should have binding energies strongly low- 
ered by shell effects, forming a region of increased shell 
stability 0,0. 

The use of "hot fusion" reactions with the neutron- 



rich 48 Ca beam and actinide targets in Dubna resulted 
in detection of 48 new nuclides with Z = 104 — 118 and 
A = 266 — 294 [7]. Several a-decay chains seen in Dubna 
were independently verified [8] . The most significant out- 
come of these recent measurements is the observed in- 
crease of half-lives with the increasing neutron number 
- consistent with the predicted increased stability of SH 
nuclei when approaching N = 184. However, the un- 
ambiguous identification of the new isotopes still poses 
a problem because their a-decay chains terminate by SF 
before reaching the known region of the nuclear chart. 
The understanding of the competition between a-decay 
and SF channels in SH nuclei is, therefore, of crucial im- 
portance for our ability to map the SH region and assess 
its extent. 

The stability of heavy and SH nuclei is profoundly af- 
fected by nuclear deformability through the competing 
fission valleys having different geometries. The optimal 
trajectory in a multidimensional space of collective co- 
ordinates that minimizes the collective action can be as- 
sociated with sequences of intrinsic symmetry-breaking 
transitions. The effects due to breaking of axial symme- 
try are known to be important around the first saddle 
[9l4l2| , and also around the second barrier [13|, LL4| . The 
reflection-asymmetric mode usually contributes at larger 
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elongations, beyond the first barrier 
intrinsic symmetry of the final system - essential for de- 
termining the final split - depends on the geometry of 
the post-saddle and pre-scission configurations of the nu- 
cleus. 

The main objective of this work is to perform realistic 



2 



predictions of decay modes of SH nuclei using an SMF ap- 
proach based on the superfluid nuclear density functional 
theory (DFT) at the deformed Hartree-Fock-Bogoliubov 
(HFB) level. The advantage of this method is its abil- 
ity to properly treat the self-consistent interplay between 
the long-ranged electrostatic repulsion and short-ranged 
nuclear attraction that gives rise to Coulomb frustration 
[4| . Our calculation - based on a realistic density depen- 
dent effective interaction between nucleons and the mi- 
croscopic description of the collective action - provides 
a quantitative description of decay properties of known 
major and minor actinides. This gives us some confidence 
in extrapolations to yet-undiscovered regions of SH nu- 
clei. While several systematic studies of fission barriers 



of SH nuclei, based on both MM 11, JjJ and SMF mod- 



els [13|, UjJ LL8| have been carried out, fission barriers are 
not observables that can be directly related to experi- 
ment. Moreover, no microscopic survey of SF properties 
on SH nuclei exists in the literature, except for some MM 
studies [Io|,[l9-22] carried out in constrained deformation 
spaces and lacking crucial self-consistent polarization ef- 
fects, and recent SMF [23|, Hif calculations limited by 
symmetry constraints imposed for most nuclei studied. 
As we shall point out in this work, imposing axial and/or 
space inversion symmetry could result in overestimation 
of SF half-lives by many orders of magnitude. 

Model — The phenomenon of fission can be understood 
in terms of many-body tunneling involving mean fields 
with different intrinsic symmetries {2^. For SH nuclei, 
the theoretical tool of choice is the self-consistent nuclear 
density functional theory (DFT) [26]. The advantage 
of DFT is that, while treating the nucleus as a many- 
body system of fermions, it provides a line for identifying 
the essential collective degrees of freedom and provides 
a starting point for time-dependent extensions [27]. To 
describe the quantum-mechanical motion under the col- 
lective barrier, it is convenient to employ the adiabatic 
time-dependent HFB (ATDHFB) theoryjp, H] that has 
been successfully applied to fission [30I, [3l|7 

The Skyrme-HFB calculations were carried out us- 
ing the framework previously discussed in Refs. 12, 3l[ 
based on the symmetry unrestricted DFT solver HFODD 
(32| capable of breaking all self-consistent symmetries 
of nuclear mean fields on the way to fission. The nu- 
clear energy density functional was approximated by the 
SkM* functional [33] in the particle-hole channel. This 
functional provides very reasonable results for fission 
barriers and SF half-lives of even-even actinide nuclei 
(34l [35| . In the particle-particle channel, we employed 
the density-dependent mixed pairing interaction [36| . To 
truncate the quasi-particle space of HFB, we adopted 
the quasiparticle-cut-off value of 60MeV in the equiva- 
lent energy spectrum [37} . As discussed in Refs. [38], such 
a large value of cut-off energy guarantees the stability of 
HFB results. The pairing strengths were adjusted to re- 
produce the neutron and proton pairing gaps in 252 Fm 



[HI; the resulting values are V n o = — 268.9MeVfm 3 and 
V p o = — 332.5 MeVfm 3 . The single-particle basis, con- 
sisting of the lowest 1140 stretched states originating 
from the 26 major oscillator shells, fully guarantees the 
stability of HFODD results jalj. All HFB states were 
taken to compute the mass tensor. 

To find the optimum trajectories in a multidimensional 
collective space, we constrain the nuclear collective coor- 
dinates associated with the multipole moments Qa/x, by 
which we explore the main degrees of freedom related to 
elongation (A/i = 20), reflection-asymmetry (A/i = 30), 
triaxiality (A/i = 22), and necking (A/i = 40). The driv- 
ing quadrupole moment Q20 is used only as a suitable 
parameter enumerating consecutive points of the one- 
dimensional collective path in a multi-dimensional con- 
figuration space. 

The microscopic ingredients needed to compute the ac- 
tion integral and penetrability are: the collective poten- 
tial energy, collective inertia (mass tensor), and collective 
ground state (g.s.) energy. To calculate the potential en- 
ergy, we subtract from the total HFB energy E tot (Q 20) 
the spurious vibrational zero-point energy ZPE(Q2o) ob- 
tained using the Gaussian overlap approximation as in 
Refs. [ioj . In this work, we use the perturbative HFB 
cranking expression for the quadrupole mass parameter 
^20,20 (Q20) [H|. The collective g.s. energy is assumed 
to be £ = 0.7ZPE(Qfo). As shown in Fig. HJ the scal- 
ing factor of 0.7 improves the agreement between exper- 
iment and theory for the SF half-lives of even-even Fm 
isotopes. Finally, the penetrability has been calculated in 
WKB according to Refs. [4l| employing action integrals 
computed along the static fission pathways. 

Results — To demonstrate that our model is capable of 
describing experimental observations, Fig. [T] displays pre- 
dicted SF half-lives for even-even Fm isotopes. This is a 
challenging case as the measured values [42[ vary within 
this isotopic chain by almost 20 decades. It is satisfying 
to see a quantitative agreement between experiment and 
theory. Similar calculations performed for the major and 
minor actinides [35| also provide good reproduction of SF 
half-lives. We wish to stress that a good agreement with 
existing data is a necessary condition for any model to 
carry out an extrapolation into the unknown region of 
SH nuclei. 

The even-even superheavy nuclei with 108 < Z < 126 
and 148 < N < 188 can be divided into three groups 
according their g.s. properties 0, IH, Q : (i) nuclei with 
prolate-deformed shapes (Q20 ~ 30 b) for N < 170; 
(ii) spherical nuclei for N > 180; and (iii) weakly de- 
formed, often triaxial systems lying between (i) and (iii). 
The nuclei with N > 180 are most stable against SF; 
they have two-humped barriers with the inner saddle at 
Q20 ~ 50 b that is higher than the outer one (Ea > 
Eb)- In most cases, triaxiality substantially reduces Ea 
0, Oil- Furthermore, for the reflections-symmetric fis- 
sion pathways with elongated fragments (sEF), triaxial- 
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FIG. 1. (Color online) SF half- lives of even-even Fm isotopes 
with 236<A<266, calculated in this study (th-0.7) compared 
with experimental data The corresponding collective 

ground state energies Eo = 0.7ZPE(Q2q) are shown in the 
lower panel. The scaling factor of 0.7 improves the agreement 
with experimental data. The results obtained without scaling 
(th-1.0) are also shown. 



ity may also reduce Eb [HI, [HI. Typically, the reflection- 
asymmetric fission valley corresponding to asymmetric 
elongated fragments (aEF) branches away from the sEF 
pathway at Q20 > 80 b beyond the inner saddle. For 
nuclei with A > 280 and Z > 108, the outer barrier van- 
ishes along aEF. SF half-lives of weakly deformed nuclei 
from the transitional region (iii) were always calculated 
relative to the prolate-deformed g.s. Both sEF and aEF 
fission valleys are taken into account in our calculations. 
The resulting fission probabilities are combined to give 
the estimated SF half-life; the larger penetrability deter- 
mines the SF mode. 

To illustrate the competition between sEF and aEF 
fission pathways, Fig. [2] shows the case of the spherical 
nucleus 306 122. The energy curves along the reflection 
symmetric fusion (sFu) and asymmetric fusion (aFu) val- 
leys are also presented. The energy gain due to triaxi- 
ality in the region of the first and second saddle can be 
assessed from the energy curves shown in the inset: the 
inner barrier is reduced by ~3 MeV by triaxiality, and 
the effect around the second saddle is weaker, around 
1 MeV. However, the outer barrier vanishes altogether 
along aEF and this favors the reflection-asymmetric fis- 
sion mode in 306 122. The total density distributions at 
pre-scission configurations in aEF and sEF are shown at 
Q20 ~ 370 b and Q20 ~ 650 b, respectively. While the 
neck rapidly vanishes in aEF, the symmetric pre-scission 
region is characterized by an extended neck. Figure [2] also 
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FIG. 2. (Color online) Total energy (circles, left scale) and 
the quadrupole mass parameter (diamonds; right scale) calcu- 
lated along sEF and aEF fission pathways in 306 122, together 
with the corresponding shapes. The energy curves along sFu 
and aFu fusion valleys are also indicated. To illustrate the 
effect of triaxiality on the inner and outer barrier, the axially 
symmetric sEF fission pathway is marked by open circles. The 
deformation energy, normalized with respect to the total g.s. 
energy, is shown in the inset. 



shows the mass parameters £?20,20 along sEF and aEF 
pathways. The two £>20, 20 trajectories are fairly similar, 
which indicates that it is the potential energy (in partic- 
ular, barrier width and height) that determines the opti- 
mal fission pathway in this case. The SF half-life along 
the axially symmetric sEF pathway is T s f — 10 13 82 s. 
Triaxial effects along sEF reduce it to 10 9,39 s, and the 
inclusion of reflection-asymmetric shapes (aEF) brings 
the SF half-life of 306 122 down to T sf = 10 6 22 s, which 
corresponds to an overall reduction of T s f by about seven 
orders of magnitude. 

The survey of the competition between sEF and aEF 
SF modes is displayed in Fig. [3l The sEF mode domi- 
nates for the Hs isotopes, SH nuclei with A < 280, and 
in a triangle defined by 290 Ds, 298 F1, and 298 Ds. For 
the remaining nuclei, the asymmetric mode is expected 
to win. In very heavy nuclei around N = 188, the bi- 
modal fission is predicted. In Fig. [3j the nuclei for which 
I log 10 (T g gp/T a j£p)| < 0.3 are marked by triangles. The 
barrier heights along aEF and sEF are similar; hence, 
it is the barrier width that determines the dominant SF 
mode. 

Figure 0] summarizes our predictions for SF and al- 
pha decay modes of SH nuclei. The calculated SF half- 
lives are shown in Fig. Ufa). The maximum value of 
T s f = 10 7 76 s corresponds to 298 F1 and T s f values of 



300 Lv, and 3 2oUdn are similar. The shortest SF half-lives, 
reaching down to 10 -10 s, are predicted for nuclei from a 
narrow corridor formed by 280 Hs, 284 F1, and 2 f|Uuo that 
lies on the border of weakly-deformed SH nuclei that ex- 
hibit prolate-oblate coexistence effects, or g.s. triaxiality 
This corridor of fission instability separates 
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FIG. 3. (Color online) Competition between sEF and aEF SF 
modes in even-even SH nuclei. The bimodal SF is expected in 
nuclei with | log 10 (T g gp/T a gp)| < 0.3 marked by coexisting 
triangles. The experimentally observed nuclei are indicated. 
The contours show the predicted SF half-lives in logarithmic 
scale: log 10 (T s //s). 
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FIG. 4. (Color online) Summary of our SkM* results for decay 
modes of SH nuclei, (a) SF half- lives log 10 T s f (in seconds), 
(b) a-decay half- lives log 10 T a (in seconds), (c) Dominant 
decay modes. If two modes compete, this is marked by coex- 
isting triangles. 



the regions of SH nuclei created in hot- and cold-fusion 
reactions. 

It is instructive to compare our SF half-lives with other 
predictions. The MM calculations [21] yield SF half-lives 
that overshoot our results by more then five orders of 
magnitude. We attribute this to the assumption of axial- 
ity and reflection-symmetry used in their work. Likewise, 



the axially symmetric HFB+D1S calculation [24| overes- 
timates our SF half-lives by many orders of magnitude. 
On the other hand, there is a nice consistency between 
our aEF results and those obtained in the axial Skyrme 
HF+BCS approach of Refs. [HI with SV-min and SV-bas 
functionals. In particular, the corridor of the maximum 
SF instability is predicted similarly by both approaches. 
It is anticipated, however, that the inclusion of triaxiality 
is likely to reduce their SF half-lives significantly. 

To estimate a-decay half-lives, we used the standard 
Viola-Seaborg expression [45] with the parameters from 
Ref. [46] and calculated Q a values. Figure SJb) sum- 
marizes our results. In general, we obtain a fairly good 
agreement with experiment. Our model underestimates 
experimental Q a values in the vicinity of the deformed 
shell closure N — 162. In this respect, the models of 
Refs. 21, lij are closer to the data. 

Our survey of lifetimes of even-even SH nuclei is sum- 
marized in Fig. HJc). According to our model, the re- 
gion of long-lived SH nuclei is expected to be centered on 
294 Ds with a predicted total half-life (considering SF and 
alpha decay) of 10 513 s, i.e.,~1.5 days. For comparison, 
the total half-life 292 Ds predicted in Refs. [UHl is 51 y 
and 14 y, respectively. 

In summary, we carried out self-consistent Skyrme- 
HFB calculations to predict main decay modes of even- 
even SH nuclei with 108 < Z < 126 and 148 < N < 188, 
assess their lifetimes, and estimate the center of enhanced 
stability in the SH region. In our model, fission pathways 
in the collective space are not constrained by imposed 
self-consistent symmetries, and ATDHFB cranking fis- 
sion inertia and zero-point energy corrections have been 
obtained microscopically. The model satisfactorily re- 
produces experimental T s f values in even-even actinides, 
which is a necessary condition for a model-based extrapo- 
lation to unknown SH nuclei. We wish to emphasize that 
our survey is the first systematic self-consistent approach 
to SF in SH nuclei that is free from artificially-imposed 
symmetry constraints that are likely to affect previous 
predictions. 

We predict two competing SF modes in SH nuclei: 
the reflection-symmetric mode sEF and the reflection- 
asymmetric mode aEF. The latter one is expected to 
prevail for N > 166, while sEF shows up in the re- 
gion of light SH nuclei and also for neutron-rich nuclei 
with N ~ 188. The region of asymmetric fission roughly 
corresponds to the region of the highest SF barriers and 
longest SF half-lives. The predicted SF half-lives of even- 
even transitional nuclei around 284 Cn are dramatically re- 
duced as compared to the current experimental estimates 
[7[. Since those systems belong to the region of shape 
coexistence and/or oblate g.s. shapes, some further in- 
crease of SF half-lives is anticipated due to the lowering 
of g.s. energy due to the shape mixing and/or appear- 
ance of a triaxial saddle at low-deformations 0, [HI, |44| . 
Other improvements of the current model include dynam- 



5 



ical treatment of penetrability by considering several col- 
lective coordinates, improved energy density functionals 
[H, and the full ATHFB inertia [3l|. Work along these 
lines is in progress. 
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